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WSTĘP 


Analiza okresowości należy do najważniejszych problemów związa- 
nych z badaniem ekonomicznych szeregów czasowych. Uwzględnienie 
wahań okresowych w modelach opisujących dynamikę zjawisk ekono- 
micznych wpływa bowiem na wzrost stopnia dopasowania tych modeli 
do danych empirycznych 1. Ponieważ poddawane dotychczas badaniu eko- 
nomiczne szeregi czasowe charakteryzowały się wahaniami o różnych 
okresach 2, szczególnego znaczenia nabiera zagadnienie identyfikacji okre- 
su wahań. Bardzo przydatną do tego celu wydaje się być analiza pro- 
cesu losowego w dziedzinie częstotliwości. Z wymienionych powodów 
w pracy omówiono najmocniejszy test okresowości prostej, czyli test 
Fishera. Podjęto również próbę zwrócenia uwagi na rolę funkcji spek- 
tralnej zarówno we wstępnej ocenie okresu wahań, jak i w ocenie istot- 


1 Zob. np. Z. Zieliński: Metody analizy dynamiki i rytmiczności zjawisk. 
gospodarczych. PWN, Warszawa 1979. 

2 W warunkach gospodarki kapitalistycznej analizowano m.in. szeregi czasowe 
z następującymi cyklami: 40—60 lat (cykl Kondratieffa), 20—30 lat (cykl Kuznetsa), 
15—20 lat (cykl charakterystyczny dla budownictwa niektórych krajów) i 2—4 lat 
(cykl Kitchina). W szeregach czasowych opisujących zjawiska ekonomiczne charak- 
terystyczne dla gospodarki socjalistycznej analizowano głównie wahania sezonowe. 
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ności wahań o różnych okresach. Dużo miejsca poświęcono problemowi 
nie omawianemu dotychczas w polskiej literaturze statystycznej, a mia- 
nowicie tzw. okresowości złożonej. Uwagi teoretyczne poparto przykła- 
dami empirycznymi wykorzystując szeregi czasowe analizowane w pra- 
cach J. Steczkowskiego, A. Zeliasia 3 i Z. Zielińskiego 4. 


KLASYCZNE METODY IDENTYFIKACJI OKRESU WAHAŃ 
WYSTĘPUJĄCYCH W EKONOMICZNYCH SZEREGACH CZASOWYCH 


W większości prac poświęconych badaniu dynamiki zjawisk ekono- 
micznych analizowane są wahania okresowe o okresie ustalonym jedy- 
nie na podstawie merytorycznej oceny zjawiska opisywanego przez dany 
szereg czasowy. Niewiele jest natomiast opracowań prezentujących me- 
tody identyfikacji okresu wahań. W polskiej literaturze statystycznej 
ciekawe uwagi odnośnie wymienionego powyżej problemu zawarto w pra- 
cach A. Sokołowskiego, К. Szymanowicz* oraz Z. Zielińskiego 6. 
W pierwszej ze wspomnianych prac dokonano identyfikacji okresu wahań 
na podstawie analizy wariancji. Omawiana metoda wymaga podziału 
п elementowego szeregu czasowego kolejno na 2, 3,...,n/2 równe części. 
Dla każdego z podziałów weryfikowana jest hipoteza o równości warto- 
ści przeciętnych odpowiadających poszczególnym podgrupom. 

Weryfikacja postawionej hipotezy dokonywaną jest zgodnie z zasada- 
mi jednoczynnikowej analizy wariancji”. Autorzy metody twierdzą, iż 
w przypadku, gdy dla danego podziału istnieją podstawy do odrzucenia 
hipotezy o równości wartości przeciętnych, można sądzić, że w szeregu 
występują wahania, o okresie równym liczbie porównywanych średnich. 

Z. Zieliński sugeruje, iż do badania okresu wahań występujących 
w ekonomicznych szeregach czasowych wystarczające jest niekiedy za- 
stosowanie testów nieparametrycznych, wśród których na uwagę zasłu- 
guje test zgodności Kendalla dla kilku zmiennych 8. 

Największe znaczenie w identyfikacji okresu wahań przypisywane jest 
testowi R. A. Fishera ®. Poniżej omówiona zostanie istota tego testu. 


3J. Steczkowski, A. Zeliaś: Analiza wariancyjna i kowariancyjna 
w badaniach ekonomicznych. PWN, Warszawa 1982. 

4 Zieliński: op. cit. 

5 А, Sokołowski, К. Szymonowicz: Analiza wariancyjna w bada- 
niach struktury hafmonicznej szeregu czasowego. Folia Oeconomica, XIX, 1976. 

s Zielinski: ор. cit. 

7 Zob. np. Steczkowski, Zelias: op. cit., s. 134—147. 

s Zob. np. Zielinski: op.cit. s. 178—180. 

® Test ten pochodzi z prac R.A. Fishera napisanych w latach 1929, 1930, 1940. 
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Załóżmy, że szereg czasowy 4x,; t=1,2,.., N} opisywany jest przy 
pomocy następującego modelu: 


x,=f(t)+u, t=1, 2, wey N (1) 
gdzie: 
f(t) — funkcja okresowa z okresami będącymi podzielnikami licz- 
by obserwacji N, 
о, — błąd losowy spełniający warunki: 


E(u) =0, Е(ы2) — 5, E(uu,) =0 dla ts 


Z kryterium Weierstrassa 10 wynika, iż funkcję f(t) można wyrazić na- 
stępująco: 


f(t) = a+)" [a, cos (21 jt/N)+Bj sin (21 jt/N)} , (2) 
gdy N=2n+1- u ‚ 
lub 
f() = a+) [а, cos (21jt/N)+B, sin (211jt/N)]-+an;(— 1)", (3) 
gdzie N=2n. = 


Amplituda drgan odpowiadajacych j/N opisywana jest w omawianym 
przypadku wzorem o postaci: 


oF +B; > j= 1, 2, өө ń (4) 
gdzie: 
R; — amplituda odpowiadająca j-tej częstotliwości, 
aj, В; — współczynniki wielomianów (2) i (3). 


Warto zaznaczyć, że harmonika z okresem 2 uwzględniona zostaje je- 
dynie w przypadku parzystej liczby obserwacji. Т. W. Anderson !! doko- 
паї oceny parametrów da, ам, оу, В; metodą najmniejszych kwadratów, 
uzyskując następujące zależności: | | 


a, =N" Уз х(, (5) 


10 Zob. А. Smoluk: Podstawy teorii aproksymacji i s.-funkcje, PWE, War- 
szawa 1974. 

п Т. W. Anderson: The Statistical Analysis of Time Series. John Wiley 
and Sons, Nowy Jork 1971. 


192 A. Góral 


N 
а, = 237 Ух, cos(2TIjt/N), j=l, 2, ..., n (6) 
t=1 
N 
b,=2N"1) x, sin(2TIjt/N), ј = 1, 2, в (1) 
t=] 
N N 
AN =N"Vx(-1y. 
t=1 


Na podstawie wzorów (4), (6), (7) uzyskiwana jest ocena amplitudy 
R; w formie: 


` 


R, = ya? +o? . j=—1,2,.. n (8) 


Hipotezę, że szereg czasowy (x,; t=1, 2,..., N} charakteryzuje się bra- 
kiem wahań okresowych można zapisać w postaci: 


H,: f(1)=£(2) ... =f(N) (9) 
lub 
Ho : R; =R= ove =R,=0. (10) 
Jezeli hipotezy alternatywne przyjmą formę: 
а (1=1, 2, ..., п; 
H,:R?>0, pozostałe В? = 0, у=) 


to mówimy o badaniu okresowości prostej 12. W przypadku, gdy w hipo- 
tezie alternatywnej zakłada się, iż amplitudy różnią się istotnie od 0 dla 
dwóch lub więcej częstotliwości, to mamy do czynienia z analizą okre- 
sowości złożonej 13. Hipotezę taką można przykładowo przedstawić w na- 
stępujący sposób: 


(j=1, 2, ..., n) 
H,:R?>0, R2>0, Л Ю?=0. (k=1,2,...,n) 
bej, k (1=1, 2, ..., п) 


Т. №. Anderson 14 wykazał, że najmocniejszym testem okresowości dla 
hipotezy głoszącej okresowość prostą jest test Fishera. Podstawę tego 
testu stanowi przy uwzględnieniu wcześniejszych oznaczeń statystyka 
o postaci: 

S=max Yj, " (11) 


1<і<п 


iż W języku angielskim: simple periodicity. 
13 W języku angielskim: compound periodicity. 
M Anderson: op. cit. 
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gdzie: 
a n a 
=R} у R. (12) 
i=l 


Hipoteza zerowa jest w omawianym przypadku odrzucana, gdy wy- 
znaczona z wzoru (11) wartość statystyki S przekroczy odczytaną z tabli- 
cy 1 wartość krytyczną gr. 

Konieczność analizy okresowości złożonej wynika z często spotykanej 
definicji wahań okresowych 15, a mianowicie: wahania okresowe o okre- 
sie P jednostek czasu stanowią sumę wahań periodycznych, których naj- 
krótszą wspólną długością cyklu jest P jednostek czasu. 

A. Siegel 16 podkreśla, że praktycznie nie ma powodów, aby sądzić, iz 
test Fishera charakteryzuje się wysoką mocą w warunkach hipotezy al- 
ternatywnej dotyczącej okresowości złożonej. Autor proponuje więc, aby 
badania nad okresowością złożoną prowadzić na podstawie statystyki 
utworzonej z wszystkich większych od gr wartości Y;. Ponieważ z teore- 
tycznego punktu. widzenia łatwo jest podać przykład, w którym występo- 
waniu okresowości złożonej towarzyszy brak według testu Fishera okre- 
sowości prostej, A. Siegel zwraca uwagę, iż rozważane zagadnienie może 
być badane na podstawie następującej statystyki: 


= У 0-0), (13) 
j=1 
gdzie: 
à — stała określona przez autora testu na podstawie ba- 
dań heurystycznych, 
(Y;—Ags)+ — oznacza max [(Y;—Agr), 0] 


Hipotezę o braku wahań okresowych odrzuca się na rzecz hipotezy 
o okresowości złożonej, gdy T przekroczy odczytaną z tabeli 1 wartość 
krytyczną. 

Dla A=1 statystyka T może być uważana za modyfikację testu Fishe- 
ra do badania okresowości złożonej. Efektem wielu badań przeprowadzo- 
nych przez E. Siegela na szeregach sztucznie generowanych jest stwier- 
dzenie, że w praktyce najlepsze rezultaty uzyskiwane są dla A=0,4. Wy- 
znaczone przez wspomnianego autora krzywe mocy testu w zależności 
od A zwracają uwagę na fakt, że gdy hipoteza alternatywna zakłada okre- 
sowość złożoną, moc testu jest znacznie wyższa dla A=0,4 niż dla A=1. 


5 Zob. S. Giembicki: Wybrane problemy analizy ekonomicznych szere- 
gów czasowych. GUS, Warszawa 1974, s. 38. 

10 A. Siegel: Testing for Periodicity in a Time Series. Journal of the Ame- 
rican Statistical Association, nr 370, 1980. 
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Tab. 1 Wartości krytyczne tp dla Tr i gy dla Sy 
The critical values tr for Ty and gp for Sy 


Pes п | er | 08er | toe | O62r | toa | Ofer | toy 
istotnosci 

5 0,684 0,547 0,137 0,410 0,274 0,274 0,412 

6 0,616 0,493 0,123 0,370 0,246 0,246 0,381 

7 0,561 0,449 0,112 0,387 0,225 0,224 0,356 

8 0,516 0,413 0,103 0,309 0,208 0,206 0,334 

9 0,477 0,382 0,0955 0,286 0,193 | 0,191 0,316 

10 0,445 0,356 0,0891 0,267 0,181 0,178 0,301 

a = 0,05 15 0,335 0,268 0,0673 0,201 0,140 0,134 0,247 

20 0,270 0,216 0,0546 | 0,162 0,116 0,108 0,213 

25 0,228 0,182 0.0462 0,137 0,0997 0,0912 0,190 

30 0,198 0,158 0,0402 0,119 0,0880 0,0791 0,173 

35 0,175 0,140 0,0357 | 0,105 0,0791 | 0,0701 | 0,160 

40 0,157 0,126 0,0322 0,0944 0,0721 0,0630 0,150 

50 0,131 0,105 0,0270 | 0,0788 | 0,0616 | 0,0525 | 0,133 


5 | 0,789 0,631 0,158 0,473 0,315 0,315 0,473 

6 | 0,722 0,577 0,144 0,433 0,289 0,289 0,433 

7 | 0,664 0,532 0,133 0,399 0,266 0,266 0,399 

8 | 0,615 0,492 0,123 0,369 0,246 0,246 0,372 

9 | 0,573 0,458 0,115 0,344 0,229 0,229 0,349 

10 | 0,536 0,429 0,107 0,322 0,214 0,214 0,329 

a = 0,01 15 | 0,407 0,326 0,0814 | 0,244 0,164 0,163 0,262 
20 0,330 0,264 0,0660 | 0,198 0.134 0,132 0,222 


25 0,278 0.223 0,0557 0,167 0,114 0,111 0,194 
30 0,241 0,193 0,0484 0,145 0,0993 0,0965 0,174 
35 0,213 7 0,171 0,0428 0,128 0,0884 0,0854 0,159 
40 0,192 0,153 0,0385 0,115 0,0799 0,0766 0,146 
50 0,160 0,128 0,0321 0,0957 0,0673 0,0638 0,128 


Żródło: A. Siegel: Testing for Periodicity in a Time Series. Journal of the Ame- 
rican Statistical Associaton, nr 370, czerwiec 1980. 


Dużą rolę w identyfikacji okresu wahań występujących w ekonomicz- 
nych szeregach czasowych może odgrywać również tzw. funkcja spek- 
tralna. Możliwości wykorzystania tej funkcji do analizy okresowości pro- 
cesów losowych zostaną omówione w następnej części pracy. | 


FUNKCJA SPEKTRALNA 
W BADANIU OKRESOWOŚCI PROSTEJ I ZŁOŻONEJ 


Niech {x,; t=1,2,...,n} oznacza realizację stacjonarnego w szerszym 
sensie i ergodycznego procesu losowego (x; t=0, + 1, +2, spoż funkcją 
autokowariancji o postaci: z 
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y(t) =E|X—M)X,,—M]|, т=0, +1, +42, ... (14) 
gdzie: 
. 1 — rząd funkcji autokowariancji, 
M — wartość oczekiwana procesu {X;}. 


Pod pojęciem funkcji spektralnej (widmowej) procesu {X,} rozu- 
miana jest transformata Fouriera funkcji autokowariancji, czyli: 


р) = У үе nik, —1/2<f<1/2 (15) 


. k=— 0 


gdzie: p(f) oznacza wartość funkcji spektralnej dla częstotliwości f. 


Zgodny i asymptotycznie nieobciążony estymator funkcji p(f) przed- 
stawiany jest często w następujący sposób: 


(fi) = 7(0)-+2 ) YC) w (k) cos(2II£jk) , (16) 
k=1 
gdzie: 
y(k) oznacza wartość estymatora funkcji autokowariancji w punk- 
cie k uzyskaną z wzoru o postaci: 
2 п-к = _ 
‚ Ү®) = У.) уж) j (17) 
=1 
- <A п, 
х=п x, (18) 
© m — punkt odcięcia funkcji autokowariancji (k= 
=0, 1, ..., m), , 
w(k) — funkcja wagowa, której transformata Fourie- 
ra nazywana jest oknem widmowym, 
f;=j/2m (j=0, 1, .., m) — częstotliwość, dla której wyznaczana jest war- 


tość widma. 


Zapewnienie wysokiej efektywności estymacji funkcji spektralnej wy- 
maga szeregu czasowego o dużej liczbie obserwacji (п > 100) oraz wła- 
ściwego doboru wartości punktu odcięcia m i funkcji wagowej w(k) 17. 

Wykorzystanie. funkcji spektralnej do analizy okresowości poprzedzo- 
ne było przez długi okres zastosowaniem do tego celu tzw. funkcji pe- 


11 Problemy te omawiane są m.in. w pracy G. M. Jenkins.a: General Con- 
siderations in the Analysis of Spectra, Technometrics, vol, 3, nr 2, maj 1961. 
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riodogramowej 18, Warto zaznaczyć, iż funkcja periodogramowa, która 
jest niezgodnym estymatorem funkcji spektralnej, może być przedsta- 
wiona w następujący sposób: 


n 2 
ы =a] X кетте —1/2<{<1[2 (19) 
t=1 
lub 
a a= 1 a 
I = 0+2 3 y(k) соѕ(2ПЕК), —1/2<f<1/2 (20) 
k=1 
gdzie: 


{x,; t=1,2,..,n) oznacza analizowany szereg czasowy, Ё jest czę- 
stotliwością, 

sy(k); k=1,2,..,n—1 oznacza ocenę wartości funkcji autokowa- 
riancji w punkcie k. 


O występowaniu wahań okresowych w badanym szeregu czasowym 
świadczą istotne wartości periodogramu dla odpowiednich częstotliwości. 

Odpowiedź na pytanie dlaczego funkcja widmowa spełnia tak istotną 
rolę w analizie okresowości procesów losowych można znaleźć m. in. 
w pracach L. Dziembały i К. Zadory 9, С. №. J. Grangera %, oraz M. Ner- 
love'a 21, 

L. Dziembała i K. Zadora podkreślają, iż analizę okresowości procesu 
losowego sprowadza się do poszukiwania na wykresie funkcji widmowej 
takich miejsc, w których gwałtownie wzrasta wartość tej funkcji. Jeżeli 
taki wzrost jest obserwowany w punktach fj, fo, ..., fk, to można wniosko- 
wać, iż rozważany proces jest cykliczny o okresach 1/fj, 1/fą, ..., L/f,. 
W przypadku, gdy w całym przedziale (0,1/2) funkcja p(f) jest bardzo 
gładka, to w badanym procesie nie występują wahania cykliczne. 

M. Nerlove 22 stwierdza, że sezonowość można zdefiniować jedynie 
przy wykorzystaniu pojęcia „funkcja spektralna”. Wspomniany autor 
pod pojęciem sezonowości rozumie tę charakterystykę szeregu czasowe- 
go, która powoduje wzrost wartości funkcji spektralnej dla częstotliwo- 


в D, R. Brillinger podkreśla w pracy pt. Wriemiennyje riady (MIR, Moskwa 
1980), że Shuster zastosował analizę periodogramową do badania okresowości już 
w latach 1894 i 1897. 

81. Dziembała, K. Zadora: Zastosowanie analizy widmowej do ba- 
dania wahań cyklicznych. Przegląd Statystyczny, nr 19, 1971. 

® С. W. J. Granger: The Typical Spectral Shape of An Economic Va- 
riable. Econometrica, vol. 34, nr 1, styczeń 1966. | 

21 M. Nerlove: Spectral Analysis of Seasonal Adjustment Procedures. Eco- 
nometrica, vol. 32, nr 3, lipiec 1964. 

2 Nerlove: op. cit. 
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ści sezonowych. С. W. J. Granger 23 omawiając problem typowej dla eko- 
nomicznych szeregów czasowych krzywej spektralnej, formułuje prawo, 
które również świadczy o tym, iż okres wahań znajduje odzwierciedlenie 
w wartościach widma. Prawo to brzmi następująco: ,,(...) jeżeli dokona- 
my dekompozycji występujących w ekonomicznych szeregach czasowych 
wahań długookresowych na składowe częstotliwościowe, to amplitudy 
tych składowych będą łagodnie zmniejszały się wraz ze zmniejszaniem 
-się okresu”. 

Warto w tym miejscu zwrócić uwagę na fakt, iż mimo że rola 
analizy częstotliwościowej w ocenie wahań okresowych jest niepodważal- 
na, to większość tego typu badań ma charakter subiektywny. Często wy- 
korzystywane jest bowiem stwierdzenie, że w szeregu czasowym wystę- 
pują wahania cykliczne w okresie T, gdy wartość funkcji spektralnej lub 
periodogramowej w punkcie 1/T powoduje istnienie w wykresie wymie- 
nionych funkcji tzw. „wierzchołka”. Wydaje się, że nadanie tego typu 
badaniom obiektywnego charakteru należy ściśle powiązać z problemem 
badania istotności wartości funkcji spektralnej. Problem ten omówiony 
został już w 1961 roku przez Е. J. Hannana 24. Wymieniony autor podjął 
próbę weryfikacji hipotezy o niezależności zmiennych tworzących dany 
proces losowy. Ponieważ krzywa spektralna jest w przypadku procesu 
czysto losowego linią prostą równoległą do osi częstotliwości, weryfiko- 
wana hipoteza jest równoważna hipotezie (10) z testu Fishera. E. J. Han- 
nan 25 wykazał, że statystyka o postaci: ` 


KE) = [ENPE] У [LEPE] 21) 
jal 
gdzie: 
fi=j/2m j=1,2,..,m 
I„(f;) — wartość funkcji periodogramowej w punkcie fi, 
p(f;) — wartość funkcji spektralnej w punkcie fj, 
charakteryzuje się rozkładem asymptotycznie zbieżnym do rozkładu sta- 


tystyki Y; z testu Fishera. Podstawę podejmowanych na podstawie roz- 
ważanego testu decyzji stanowią wartości następującej statystyki: 


Ѕ = max K(,) = max л ГУ [LEDE] | (22) 


j=l 


gdzie: р(ї)) oznacza ocenę wartości funkcji spektralnej w punkcie fj. 


3 Granger, op. cit. 

aE. J. Hannan: Testing for a Jump in the Spectral Function. Journal of 
the Royal Statistical Society, B, vol. 23, nr 2, 1961. 

8 Е J. Hannan: Analiz wriemiennych riadow. Nauka, Moskwa 1964. 
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W przypadku gdy uzyskana na podstawie określonego szeregu czaso- 
wego wartość statystyki Sp przekracza odczytaną z tablicy 1 wartość 
krytyczną gy, hipotezę o braku wahań okresowych należy odrzucić. 

Przytoczone w drugiej części pracy twierdzenie T. W. Andersona od- 
nośnie mocy testu Fishera dowodzi, iż test ten powinien odgrywać de- 
cydującą rolę w badaniu okresowości prostej. Wydaje się natomiast, że 
związana z funkcją spektralną statystyka K(f;) może być pomocna przy 
konstrukcji testu służącego do weryfikacji hipotezy o okresowości złożo- 
nej. 
Wykorzystanie uwag A. Siegela odnośnie statystyki Ty prowadzi do 
wniosku, iż okresowość złożoną można badać na podstawie statystyki 
wyrażonej w następujący sposób: 


Т; = >. |К()—Лвь]+ » (23) 
j=l 


gdzie: A zgodnie z wcześniejszymi uwagami przyjmuje się na poziomie 
0,4. 


Ponieważ rozkład T, jest asymptotycznie zbieżny do rozkładu Ty, hi- 
poteza o braku wahań okresowych w danym szeregu czasowym jest od- 
rzucana, gdy wartość statystyki Т„ przekracza odczytaną z tablicy 1 war- 
tość krytyczną. 


BADANIA EMPIRYCZNE 


Przedstawiony na ryc. 1 szereg odchyleń od trendu wskaźników wy- 
korzystania taboru PKS w Szczecinie w poszczególnych miesiącach lat 
1965—1972 analizowany był przez Z. Zielińskiego **. Autor ten weryfi- 
kował hipotezę o występowaniu w wymienionym szeregu wahań sezo- 
nowych na podstawie nieparametrycznego testu Kendalla. Test ten po- 
zwolił stwierdzić istotność analizowanych wahań. Analogiczne badanie 
przeprowadzone zostanie przy wykorzystaniu testów: Fishera i Hannana. 
Postawiono następujące hipotezy: 


H, M R,;=Re= aes Rn, n=48 


H,:R;>0, Ri=0 A j=l, 2,..., п. 
iźj 
Na podstawie wzorów: (5), (6), (7) i (9) wyznaczono dla każdego j= 
=l,2,...n oceny amplitud Rj. Wartość statystyki S, która jest spraw- 


26 Zielinski: op. cit., s. 180. · 
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Ryc. 1. Odchylenie od trendu wskaźników wykorzystania taboru PKS w Szczecinie 
w latach 1965—1972 
Deviations from trend of coefficients of using Polish Motor Transport 
| (PKS) in Szczecin in 1965—1972 


dzianem w teście Fishera, uzyskano z zależności (11) i (12) przyjmując 
n=48. Wartość ta wyniosła 0,27. Ponieważ S przekracza odczytaną z ta- 
blicy 1 dla a=0,05 i n=48 wartość krytyczną gy (gr 0,14), to hipote- 
ze Н, należy odrzucić na rzecz alternatywy Ну, gdzie 77 j=8. Widać więc, 
że test Fishera potwierdza sąd Z. Zielińskiego, iż w szeregu występują 
wahania o okresie rocznym. Prostą okresowość przedstawionego па ryc. 1 
szeregu, analizowano również przy pomocy spektralnego testu Hannana. 
Wartość statystyki Sy okazała się jednak nieistotna 28 na poziomie istot- 
ności a=0,01 i dla m=24. 

Na ryc. 2 przedstawiono szereg czasowy zawierający informacje od- 
nośnie skupu mleka przez spółdzielnie mleczarskie w woj. krakowskim 
w latach 1961—1968. J. Steczkowski і A. Zeliaś *9, wykorzystując anali- 
zę wariancji stwierdzili istotność występujących w danym szeregu wa- 
hań sezonowych i trendu. Wykonując, analogiczne jak w poprzednim 
przypadku, obliczenia uzyskano następujące wartości Sp i Sp: Sp = 0,38, 
Sy = 0,9. Wartości te przekraczają odczytane z tablicy 1 wartości kry- 
tyczne 30 gp, co potwierdza wniosek J. Steczkowskiego i A. Zeliasia. 


27 Okazało się, że S = Ү,. 

28 Sg=0,1Lra gr = 0,23. 

2 Steczkowskt Zeliaś: op. cit. 

30 W tym przypadku wartości krytyczne dla Sp oraz Sp odczytywane są odpo- 
wiednio przy założeniu, że а = 0,05, п = 48 і a = 0,05, п = 24. 
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Ryc. 2. Skup mleka przez spółdzielnie mleczarskie w woj. krakowskim w latach 
: 1961—1968 
Purchasing milk by co-operative creamery in Kraków province in 1961— 
1968 


Podobne badanie zdecydowano sie przeprowadzić dla szeregu uzyska- 
nego w wyniku eliminacji trendu występującego w tzw. szeregach jedno- 
imiennych czasokresów 31. Do eliminacji trendu wykorzystano różnicowa- 
nie poszczególnych wartości szeregu wg następującej zależności: 


к= Xi+1,k Жүн. і = 1, 2, A 7 (lata) 
j = 1, 2, ..., 12 (miesiące) 


Statystyki Sp, Sp przyjęły odpowiednio wartości: 0,11 i 0,56. Widać 
więc, że jedynie test Hannana pozwolił stwierdzić, iż w szeregu {Zik} Wy- 
stępują istotne wahania okresowe. 

Warto zaznaczyć, że w rozważanych przypadkach nie było potrzeby 
wykorzystywania testów dotyczących badania okresowości złożonej. Sto- 
sówanie tych testów zaleca się wówczas, gdy w szeregu stwierdza się 
brak wahań określonego typu. Może się bowiem okazać, że wartości sta- 


31 Zob. A. Zeliaś: Teoria prognozy. PWE, Warszawa 1979. 
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tystyk Y;, K(f;), pomimo iż są nieistotne z punktu widzenia testów Fi- 
shera i Hannana, to łącznie dają wartość Ty oraz Ts przekraczające od- 
czytaną z tablicy 1 wartość krytyczną t, 


ZAKOŃCZENIE 


O dużej roli, jaką można przypisać analizie procesów losowych 
w dziedzinie częstotliwości decyduje niewątpliwie fakt, iż wahania okre- 
'sowe o okresie Р należy traktować jako sumę wahań periodycznych 
o okresach, których najkrótszą wspólną długością cyklu jest P jednostek 
czasu. Prezentowane w pracy uwagi teoretyczne oraz badania empirycz- 
ne upoważniają do sformułowania następujących wniosków: 

1) badanie okresowości prostej należy prowadzić w oparciu o test 
Fishera, 

2) wykres funkcji spektralnej można jedynie wykorzystać do wstęp- 
nej oceny okresu wahań występujących w danym szeregu czasowym, 

3) obiektywna ocena istotności  „wierzchołków” występujących 
w krzywej spektralnej wymaga zastosowania testu Hannana, 

4) w sytuacji, gdy testy okresowości prostej prowadzą do odrzucenia 
hipotezy o występowaniu w szeregu wahań okresowych, badanie należy 
uzupełnić o analizę okresowości złożonej na podstawie testu A. Siegela 
lub testu zaproponowanego w pracy. 

Na zakończenie warto zaznaczyć, że obliczenia do niniejszej pracy wy- 
konano w ŚCO CYFRONET — Kraków na podstawie własnych progra- 
mów. 


ы РЕЗЮМЕ 


В статье обращается внимание на роль частотного анализа в исследовании 
периодичности случайных процессов. Обсуждается наиболее мощный критерий 
простой периодичности — критерий Фишера, демонстрируются достоинства 
спектральной функции в идентификации периода колебаний, наблюдаемых 
в экономических временных рядах. Большое внимание уделено He обсуждае- 
мой в польской, статистической литературе проблеме исследования сложной 
периодичности. Предлагается критерий, использующий значение спектральной 
функции для анализа этой проблемы. 

Теоретические замечания проверяются на основании анализа периодично- 
сти временных рядов, описывающих экономические явления. 


SUMMARY 
The work directs attention to the role of frequency analysis in investigating 


periodicity of stochastic processes. The most powerful test of simple periodicity has 
been discussed, that is Fisher’s test and advantages of spectral function have been 
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presented in identifying the period of fluctuations occurring in economic time 
series Much space has been devoted to. the problem to investigating complex 
periodicity, the problem which has not been discussed in Polish statistical litera- 
ture. The author proposed a test taking advantage of the values of spectral function 
for the analysis of the problem. ` 

Theoretical remarks have been verified on the basis of the analysis of perio- 
dicity of time series describing economic phenomena. 


